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Abstract We intercompare the above-cloud aerosol optical depth (ACAOD) of biomass burning plumes 
retrieved from A-train sensors, i.e., Moderate Resolution Imaging Spectroradiometer (MODIS), Cloud-Aerosol 
Lidar with Orthogonal Polarization (CALIOP), Polarization and Directionality of Earth Reflectances (POLDER), and 
Ozone Monitoring Instrument (OMI). These sensors have shown independent capabilities to retrieve aerosol 
loading above marine boundary layer clouds — a kind of situation often found over the southeast Atlantic Ocean 
during dry burning season. A systematic comparison reveals that all passive sensors and CALIOP-based 
research methods derive comparable ACAOD with differences mostly within 0.2 over homogeneous 
cloud fields. The 532 nm ACAOD retrieved by CALIOP operational algorithm is underestimated. The 
retrieved 1 064 nm AOD however shows closer agreement with passive sensors. Given the different types 
of measurements processed with different algorithms, the reported close agreement between them is 
encouraging. Due to unavailability of direct measurements above cloud, the validation of satellite-based 
ACAOD remains an open challenge. The intersatellite comparison however can be useful for the relative 
evaluation and consistency check. 


1. Introduction 

Atmospheric aerosols in cloud-free scenes generally produce a net negative radiative effect (cooling) at the 
top-of-atmosphere (TOA). However, elevated layers of absorbing aerosols above cloud can potentially exert a 
significant level of atmospheric absorption and produce a positive radiative forcing at TOA [Keil and 
Haywood, 2003; Chand et at., 2009]. The magnitudes of the forcing directly depend on the microphysical and 
optical properties of the aerosol layer and underlying cloud. Conventional satellite retrievals of aerosols are 
restricted to the cloud-free regions, which severely limits our understanding of the aerosol effects on cloud 
radiative forcing and microphysical properties. 

There has been a growing interest in the recent years on quantifying the aerosol loading above cloud from 
satellite-based active as well as passive measurements. Several algorithms have been introduced by different 
research groups employing different types of satellite measurements to retrieve the above-cloud aerosol 
(ACA) optical depth (ACAOD). The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) sensor employs 
an operational aerosol algorithm that uses estimated vertical extinction profile derived from the backscatter 
measurements and assumed lidar ratio (extinction to backscatter) to derive the AOD in clear as well as cloudy 
skies [Winker et at., 2009; Young and Vaughan, 2009]. Alternative CALIOP-based research methods, namely, 
the depolarization ratio (DR) [Hu et at., 2007] and the color ratio (CR) [Chand et at., 2008], respectively, have 
also been introduced to retrieve ACAOD. Waquet et at. [2009] make use of polarized radiances measured by 
the Polarization and Directionality of Earth Reflectances (POLDER) to retrieve AOD above cloud. An opera- 
tional version of this algorithm has also been developed, which is capable of retrieving ACAOD globally 
[Waquet et at., 2013]. Taking the advantage of enhanced aerosol absorption sensitivity in the near-UV, Torres 
et at. [2012] introduced an innovative technique to retrieve ACAOD and underlying aerosol-corrected cloud 
optical depth (COD), simultaneously, from observations by the Ozone Monitoring Instrument (OMI). Jethva 
et at. [2013] have developed a color ratio method which employs one visible and one shortwave IR (SWIR) 
channel measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS) to deduce a pair 
of ACAOD and aerosol-corrected COD, simultaneously. 
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Validation of the above-cloud AOD retrieval is a challenging task primarily due to lack of adequate direct 
measurements of aerosols above cloud. A few field campaigns such as Southern African Fire-Atmosphere 
Research Initiative 2000 [Haywood et a!., 2004] and Tropospheric Aerosol Radiative Forcing Observation 
Experiment carried out aerosol measurements in the clear as well as cloudy atmosphere. Flowever, these data 
sets have limited samplings in time and space and also not meant to measure the aerosols above cloud in 
particular. In this situation, an intercomparative analysis of the existing research retrievals from different 
sensors would be helpful to check the consistency (or lack thereof) among entirely independent techniques. 
In this paper, we intercompare the coincident research retrievals of ACAOD retrieved using six entirely in- 
dependent techniques applied to four sensors, i.e., MODIS, CALIOP, POLDER, and OMI-flying in line formation 
in the A-train satellite constellation. Two events of biomass burning aerosols above cloud are selected over 
the southeast Atlantic Ocean (SEAO), where smoke aerosols generated from the intense biomass burning in 
central Africa are often transported westward over the oceanic low-level stratocumulus cloud deck. We 
emphasize here that the present study does not constitute the "validation" of the ACA retrieval, instead, the 
objective of this study is to check the consistency among entirely independent techniques that have been 
developed to measure the above-cloud aerosol loading. A brief description of the different ACA techniques 
and data set are given in section 2. The results of intercomparison are presented in section 3. Possible sources 
of uncertainties associated with each method are discussed in section 4. 

2. A Brief Description of ACA Techniques and Data Sets 

2.1. MODIS 

The presence of absorbing aerosol layer above cloud reduces theTOA reflectance and CR between the visible 
and SWIR wavelengths. The general CR technique developed by Jethva et al. [2013] exploits this signal and 
uses reflectance at two channels, i.e., 470 nm and 860 nm to retrieve ACAOD and underlying COD, simulta- 
neously. The method requires MODIS TOA reflectance (MYD021 KM), geolocation data (MYD03), and MODIS 
cloud product (MYD06) for cloud screening all of them at 1 km resolution. These products were obtained 
from http://ladsweb.nascom.nasa.gov/data/. The aerosol model required to generate look-up-table (LUT) was 
based on the monthly statistics of Aerosol Robotic Network (AERONET) measurements at an inland site, 
Mongu (1 5°S, 23°E), Zambia, in southern Africa. The retrieved ACAOD at 860 nm was converted to 500 nm 
following the spectral extinction assumed in the aerosol model. 

2.2. OMI 

The near-UV technique applied to the OMI observations makes use of two-channel measurements (354 and 
388 nm) to retrieve ACAOD and underneath COD, simultaneously. An increased sensitivity of TOA signal to the 
aerosol absorption above cloud in the near-UV region has been exploited to retrieve ACA, which, under a pre- 
scribed set of assumption, can be directly related to the magnitude of ACAOD and COD. The OMI TOA orbital 
data for the two case studies was obtained from NASA Goddard Earth Sciences Data and Information Sen/ices 
Center (GES-DISC) (http://disc.sci.gsfc.nasa.gov/Aura/data-holdings/OMI/omaeruv_v003.shtml). Aerosol micro- 
physical properties assumed for LUT were acquired from the existing OMAERUV standard carbonaceous aerosol 
models. A full description of the algorithm is given by Tones et at. [201 2]. The retrieved near-UV ACAOD at 388 nm 
was converted to 500 nm, assuming the model-based spectral extinction for the intersatellite comparison. 

2.3. CALIOP 

The operational CALIOP algorithm employs feature and layer detection scheme [Winker et al., 2009] com- 
bined with an extinction retrieval algorithm [Young and Vaughan, 2009] that requires an assumption on the 
extinction-to-backscatter ratio to retrieve the height-resolved extinction profile and columnar AOD at 532 nm 
and 1064nm wavelengths at 5 km scale (level 2). The level 2 CALIOP data (V3.01) were obtained from NASA 
Langley's Atmospheric Science Data Center (ASDC) web portal (https://eosweb.larc.nasa.gov/project/calipso/ 
cal_lid_l2_05kmalay-prov-v3-01_table). 

In addition to the operational CALIOP method, two alternative active remote sensing techniques, namely, 
depolarization ratio and color ratio, have been developed to deduce ACAOD from CALIOP observations. 
These techniques are based on light transmission methods that use bright and commonly available lidar 
targets of liquid water clouds underneath the aerosol layers. These methods utilize three properties, i.e., layer 
integrated attenuated backscatter and depolarization ratio at 532 nm, and layer-integrated attenuated color 
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Table 1 . Salient Properties of the ACA Retrieval Techniques and Associated Sensors 

Physical Basis Algorithmic Assumptions Input Parameters Retrieved Parameters Retrieval Uncertainty 


MODIS Jethva et al. 
[2013] 


CALIOP Operational 
Method Winker et al. 
[2009]; Young 
and Vaughan [2009] 

CALIOP Depolarization 
Method Hu et al. 
[2007] 

CALIOP Color Ratio 
Method Chand et al. 
[2008] 

POLDER Waquet et al. 
[2009, 2013] 


OMI Torres et al. [2012] 


Change in the 
reflectance and color ratio 
in the visible 
(VIS)/SWIR when 
absorbing aerosols 
overlay cloud 
Vertically resolved 
aerosol backscatter 
and inferred 
extinction profile 

Change in the 
transmittance (DR) 
VIS/IR when absorbing 
aerosols overlay cloud 
Change in the VIS/IR 
transmittance 
when absorbing aerosols 
overlay cloud 
Creation of polarization at 
forward scattering angles. 
Reduction of the polarized 
signal in the cloudbow. 


Change in the near-UV 
reflectance and UV 
aerosol index (UVAI) 
when absorbing aerosols 
overlay cloud 


Aerosol/cloud size 
distribution, aerosol 
spectral extinction 
and absorption, aerosol/ 
cloud profiles, and surface 
albedo 

Extinction to backscatter 
ratio and feature/layer 
detection scheme 


Cloud extinction-to- 
backscatter ratio at 532 nm 
and Rayleigh correction 
above cloud 
Angstrom Exponent 
(532-1 064 nm) and color 
ratio (VIS/IR) of underlying 
cloud layer 

Six fine-mode spherical 
aerosol models with 
refractive index of 1.47-0.01 i, 
one nonspherical- 
mineral dust model, 
and only one aerosol/cloud 
profile. 

Aerosol/cloud size 
distribution, aerosol spectral 
extinction and absorption, 
aerosol/cloud profiles, 
and surface albedo 


TOA reflectance at 470 
and 860 nm at 
1 x 1 km 2 
resolution 


Attenuated backscatter 
profile at 532 and 1 
064 nm at 0.3 km 
scale 

Layer integrated 
attenuated backscatter 
and depolarization 
ratio at 532 nm 
Layer integrated 
attenuated backscatter 
and color ratio 
at 532 nm 

TOA polarized radiance 
at 670 and 865 nm at 
6x6 km 2 

resolution (POLDER) and 
cloud droplets effective 
radius (MODIS) 


TOA reflectance at 
388 nm 

and measured UVAI 
at 1 3 x 24 km 2 resolution 


ACAOD and aerosol- 
corrected COD at 
860 nm; converted 
to 500 nm based on 
the model assumed 
spectral extinction 
Vertical extinction 
profile, AOD at 532 
and 1064nm and at 
0.3 and 5 km (level 2) 
scales 

ACAOD at 532 nm 


ACAOD at 532 nm 
and Angstrom exponent in 
the 532-1064 nm range 

AOD at 865 nm and 
Angstrom exponent. 
The AOD is converted at 
500 nm using the retrieved 
Angstrom exponent. 


ACAOD and aerosol- 
corrected COD at 388 nm; 
converted to 500 nm based 
on model assumed spectral 
extinction 


Depends on ACAOD and 
COD. Typically between —12 
and 46% at COD of 10 and 
AOD of 0.5 for an uncertain 
single-scattering albedo 
(SSA) of ±0.03 
From Winker et al. [2009]: 
0.05 + 0.40 x AOD, 


Depends on the integrated 
attenuated backscatter at 
532 and its depolarization 
ratio. 

Depends on the integrated 
attenuated backscatter 
at 532 nm and 1064nm and 
Angstrom exponent 
Depends on AOD and 
microphysics. For an AOD 
of 0.2 (at 865 nm) : AOD error 
of 0.05 for a real refractive 
index uncertainty of +/— 0.06 
and error of 0.02 for an 
imaginary refractive index 
uncertainty of +/— 0.01 
Depends on ACAOD and 
COD. Typically between —23 
and 43% at COD of 10 and 
AOD of 0.5 for an uncertain 
SSA of ±0.03 


ratio. Both the transmission methods require some screening criteria and independent calibration constants. 
The details of the DR and CR methods are given by Hu et al. [2007] and Chand et al. [2008], respectively. 

2.4. POLDER 

Aerosols above clouds generate an additional polarized signal at forward scattering angles and reduce the 
polarized signal of the cloudbow. The operational algorithm developed for the POLDER instrument takes 
advantage of these effects to retrieve the AOD above clouds and an aerosol model. The method consists of a 
comparison between the polarized radiance measured at 670 and 865 nm and precomputed polarized 
radiances using six spherical fine-mode and one coarse mode mineral dust aerosol models. The aerosol 
properties are first retrieved at a spatial resolution of 6 x 6 km 2 (i.e., POLDER native resolution) and then 
aggregated at a resolution of 1 8 x 1 8 km 2 . The cloud properties retrieved from MODIS and POLDER 
were colocated to characterize the cloudy scenes. The algorithm is applied only to homogeneous cloudy 
pixels associated with optically thick liquid water clouds (COD > 3). More details about the various filters 
used to obtain a quality-assessed product and aerosol models can be found in the work of Waquet et al. 
[2013]. Table 1 summarizes the basic properties of different ACA techniques that are applied to the 
A-train sensors. 

2.5. Colocation 

Different sensors on the A-train satellites observe the same part of the Earth within a few minutes. The 
relative position of the A-train satellites (hence time difference) has gone through some changes since 
the inception of the constellation (http://atrain.nasa.gov/historical_graphics.php). The case studies 
examined in this paper belong to August 2006 and 2007, during which the time differences between Aqua and 
CALIPSO, CALIPSO and PARASOL, and PARASOL and Aura were about 1, 2, and 5 min, respectively. It is assumed 
that the cloud and aerosol systems do not change significantly during such short span. The spatial resolution 
of the A-train sensors however is different and therefore intercepts different coverage. While high-resolution 
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Figure 1. An above-cloud aerosol event observed over the southeast Atlantic Ocean on 13 August 2006. (top) The spatial plots of ACAOD 
retrieved from passive sensors are displayed. The CALIOP track is imposed on the MODIS plot, (bottom) A comparative plot of ACAOD along 
the CALIOP track. 


sensors such as MODIS and CALIOP make measurements at subkilometer scales, POLDER on PARASOL platform 
intercept Earth with a footprint of size 6x6 km 2 and OMI on Aura has even larger footprint of 1 3 x 24 km 2 . We 
consider the CALIOP along-track level 2 geolocation as a reference track for the intercomparison of ACAOD. In 
order to match with the 5 km averaged CALIOP retrievals, MODIS's 1 km ACAOD retrievals were also averaged to 
5 km 2 grid box. Due to their coarser spatial resolutions, OMI and POLDER ACAOD retrievals were not averaged. 
We adopt the closest pixel approach (maximum distance ~5 km) to colocate passive sensors with the lidar 
retrievals. 


3. Results 

We have selected SEAO as a test bed to intercompare the different satellite retrievals. During dry season (July 
through September), the region-wide intense biomass burning over southern Africa emits huge amounts of 
carbonaceous particulate matter into the atmosphere, which, under the influence of prevailing westerly 
winds, is often advected over the boundary layer marine stratocumulus clouds. Two case study events of 
ACAOD are selected for the present analysis. 

Figure 1 (top left) shows a true color RGB image made from two MODIS LI 5 min granules over SEAO on 13 
August 2006. The rest three plots on top display the spatial distribution of ACAOD retrieved from Aura/ 
MODIS, Aqua/OMI, and PARASOL/POLDER, respectively. All three passive sensors, in general, retrieve similar 
spatial pattern of ACAOD with maximum aerosol concentration off the western coast of southern Africa. 
Spatial coverage of ACAOD retrieved by three sensors is different. MODIS provides maximum retrieval 
availability followed by OMI and POLDER. Due to their coarser spatial resolution, OMI and POLDER algorithms 
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Figure 2. Same as in Figure 1 but for an above-cloud aerosol event observed by the A-train sensors on 2 August 2007. 


reject pixels with partial cloud cover, whereas the finer pixel size of MODIS (250-500 m) can adequately 
capture small-scale clouds, thus yield more retrieval coverage. 

Figure 1 (bottom) shows a comparative plot of ACAOD along the CALIOP track for the same case. All three 
passive sensors (MODIS, POLDER, and OMI) are found to be in close agreement (AACAOD-0.1) between the 
latitude coverage 23°S to 15°S. Between 15°S and 12°S, MODIS tends to retrieve higher AOD than those of 
OMI and POLDER. Throughout the spatial extent of this event (30°S to 5°S), CALIOP's standard 532 nm 
retrievals are lower by a factor as large as 5 or even more relative to the passive sensors. However, CALIOP 
standard 1064nm retrieval converted to 500 nm using POLDER's Angstrom exponent brings lidar retrieval 
closer to the passive sensors. 

Largest discrepancy between all sensors is seen between the equator and 5°N, where the true color image 
shows presence of small-scale broken clouds. Over this area, MODIS and CALIOP-based DR method retrieve 
very large values of AOD (—2.0) followed by POLDER and CALIOP 1064nm retrieving AOD of about 1.3; 
CALIOP-based CR method derives a range of values between 0.5 and 1, and finally CALIOP standard 532 nm 
AOD (-0.3-0.4) are falling the lowest in the comparison. A few retrievals from CALIOP and MODIS around 
latitude 6°S correspond to thin clouds. POLDER and OMI do not attempt retrieval over this region due to 
partial cloud coverage in their bigger footprints. However, owing to have higher resolution, both CALIOP and 
MODIS are able to detect the presence of aerosols above cloud over small-scale cloud fields. The magnitudes 
of retrieved ACAOD are however significantly different. 

The results for the second case study, 2 August 2007, are shown in Figure 2. Similar to the previous case, all 
three passive sensors retrieve consistent spatial distribution of ACAOD with variable spatial coverage. The 
ACAOD intercomparison plot in Figure 2 (bottom) shows close agreement (AACAOD-0.1 -0.2)) among all 
the sensors except CALIOP 532 nm retrieval between latitude 23°S and 9°S. A few retrievals from CALIOP 
and MODIS around latitude 5°S correspond to very thin and shallow clouds, where both sensors retrieve 
significantly different AOD. Also, the discrepancy between all sensors is large over the broken cloud region 
between the equator and 5°N. 
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4. Discussion and Conclusion 

The accuracy of the ACAOD retrieval from each sensor depends on several factors, which are briefly discussed 
here. All methods compared here assume that the sensor footprint is fully overcast. This requirement can be 
assured by the higher spatial resolution sensor MODIS and CALIOP. Coarser resolution sensors OMI and 
POLDER intercept larger area at surface and therefore require robust cloud mask. The impact of sensor 
resolution is reflected in the different spatial coverage of retrieval in Figures 1 and 2. Moreover, the air mass 
and thus AOD captured by different sensors directly depends on the pixel size and (in)homogeneity of 
aerosol loading. This can be a source of disagreement between different sensors. The particle microphysics 
and cloud droplet distribution assumed in the aerosol model can be an important source of uncertainty in the 
ACAOD retrieval. For instance, the color ratio algorithm applied to MODIS and OMI requires assuming the 
aerosol and cloud size distribution, and more importantly the imaginary part of the refractive index or 
single-scattering albedo (to). The aerosol models used in the CR algorithms are essentially derived from the 
AERONET inversions, which themselves carry retrieval uncertainty; for instance, co is expected to be 
accurate within ±0.03. An uncertain co beyond 0.03 can result in an ACAOD retrieval error of about 50% 
or larger [Torres et at., 2012; Jethva et al., 2013]. 

The retrieval accuracy for the operational POLDER algorithm mainly depends on the assumptions made for 
the aerosol microphysical properties. Polarization measurements at 670 and 865 nm are primarily sensitive 
to the optical thickness of scattering. The method retrieves the total (scattering and absorbing) AOD by 
assuming a constant value of 0.01 for the imaginary part of the complex refractive index. A relative AOD error 
of 1 0% at 865 nm is expected for an imaginary refractive index uncertainty of 0.01 . The assumption made for 
the real refractive index also impacts the AOD retrievals (see Table 1). The good quantitative agreement 
observed between the POLDER, OMI, and MODIS AODs, estimated at 500 nm, suggests that the Angstrom 
exponent retrieved by POLDER is also well estimated. The POLDER algorithm is only applied to homogeneous 
cloudy pixels associated with a cloud fraction of 1 00% and a COD larger than 3, thus ensuring the accuracy of 
the retrieved AOD. 

The CALIOP-based CR method assumes the aerosol Angstrom exponent and the calibration constant or color ratio 
for unobstructed cloud targets. For instance, an error of —0.5 in Angstrom exponent can cause about 18% under- 
estimation in CR-based ACAOD relative to that of DR result. Actual variations in these parameters, if not accounted 
for in the inversion, are the possible sources of uncertainty. The accuracy of DR method relies on the accuracy of the 
lidar calibration and assumed layer-integrated attenuated backscatter (532 nm) for unobstructed clouds. 

While the CALIOP standard 532 nm AOD is found to be significantly lower by a factor of 4 to 6 relative to 
passive retrievals, the 1 064 nm ACAOD when converted to the 500 nm using the POLDER-retrieved Angstrom 
exponent (black triangles in Figures 1 and 2) shows a closer agreement with retrievals from passive sensors. 
An analogous result has also been observed by Omar et al. [2013], where CALIOP standard 532 nm AOD 
retrievals were found to underestimate AERONET measurements by about 40-50%, while its 1064nm 
retrievals did not show any obvious underestimation and were in better agreement with those of AERONET. 
Kacenelenbogen etal. [201 1] also noted that CALIOP is underpredicting cloud-free AOD (532 nm) of a biomass 
burning aerosol plume by a factor of 2 compared to the coincident high spectral resolution lidar, ground- 
based, and other satellite measurements for a case study over the mid-Atlantic region in the U.S. For an 
optically thick absorbing aerosol layer, the two-way lidar signal undergoes a strong attenuation which can 
lead to weak or even loss of signal particularly from the lower layers. A possible explanation of the reduced 
signal in the presence of carbonaceous particles is the effect of aerosol absorption that, at the typically large 
optical depth values of smoke layers, would yield large aerosol absorption optical depths that significantly 
reduce the number of backscattered photons [Torres et al., 2013]. This effect would be more pronounced at 
shorter wavelengths (owing to the stronger extinction: scattering and absorption) than at longer wave- 
lengths (1064 nm). As a result, the CALIOP 532 nm signal detects only a fraction of the absorbing layers, 
whereas the 1064nm measurements observe the entire aerosol column above an opaque surface (refer to 
the curtain images of CALIOP backscatter for the present two case studies at http://www.calipso.larc.nasa. 
gov/products/lidar/browse_images/). Furthermore, misclassification of aerosol types and subsequent use of 
inappropriate lidar ratio can be an additional source of uncertainty in the CALIOP standard retrieval. 

A detailed analysis on the attribution of the differences to various sources of uncertainties for the present 
case studies is currently out of the scope of this paper due to space limitation. The sensitivity analysis of the 
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ACAOD retrievals is reported in the lead algorithm papers referenced in the Introduction. Due to the scarcity 
of direct measurements of ACAOD, it is hard to quantify the actual uncertainty in each of the present retrieval 
methods. The suborbital direct measurements of aerosols above cloud are strongly needed over the regions 
of SEAO (smoke above cloud), tropical Atlantic Ocean (smoke and dust above cloud), and Southeast Asia 
(smoke and pollution above cloud), where the elevated layers of absorbing aerosols over cloud are often 
observed from satellites. 

Given that each ACAOD method is designed independently and relies on different types of measurements 
from different sensors, an overall close agreement between them over the homogeneous cloud fields is an 
encouraging result. Such intercomparison should be further extended to multiple events, regions, and spa- 
tial/temporal scales for the comprehensive evaluation. It is expected that the ACAOD research methods will 
be applied to produce operational products in the coming years, which in conjunction with the standard 
cloud-free aerosol products will provide us with an unprecedented all-sky aerosol distribution from space. 
This can substantially enhance our knowledge on how aerosols affect the cloud radiative forcing and mi- 
crophysical properties, and aerosol transport. 
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